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APPUCATION  OF  CONTINUOUS,  ORTHONORMAL 
WAVELETS  TO  IMAGE  DATA  REDUCTION 


GENERAL  INTRODUCTION 

It  is  often  possiUe  to  adequatdy  rq^esoit  a  digital  image  usiqg  significantly  less 
information  than  is  required  to  qpedfy  each  individual  picture  element.  Sudi  a  rqnresentatitm  is 
usefiil  since  the  image  can  then  be  stored,  tranamtted,  and  di^li^ed  more  effidently  (see,  e.g., 
Geri,  Zeevi,  &  Porat,  1990;  Gertner  &  Geri,  1994).  The  dassical  tedmiques  of  Fourier  analysis 
(FA)  have  been  used  to  effidently  rqnvsent  inqages  by  a  set  of  spectral  fimctkms  whidi,  vtdmi 
added  together,  will  adequatdy  rqwoduce  the  original  image.  Fourier  analysis  howeva*,  is  not 
ideal  for  rq>resenting  natural  images  «nce  the  latter  are  generally  nonstationary  ->  that  is,  thdr 
characteristics  vary  as  a  function  of  location  within  the  image.  The  q>ectral  fiinctions  of  FA  etu^h 
encompass  the  entire  image,  and  hence  they  are  not  wdl  suited  for  rqvesenting  the  localized 
features  of  natural  images. 

Wavelet  decompo^on  involves  finding  image  expanaon  coeffidents  with  respect  to  a 
basis  derived  fi-om  those  wavelets.  The  wavelet  coeffidoits,  or  some  subset  of  them,  can  then  be 
stored,  transmitted,  compressed,  and  used  for  image  generation.  The  effident  rq)resaitation  of 
images,  thovfore,  requires  that  coeffidoits  be  discarded  such  that  the  image  generated  fi-om  the 
remaining  coeffidents  is  perceptually  equivalent  to  the  oripiul.  We  have  decomposed  images  and 
recmistructed  reduced  veraons  of  them  by  selectively  discarding  hi^-fi-equency  conq)onents  that 
are  kimwn  to  be  less  irnportam  in  conveying  visual  form  information.  The  reduced,  fidl 
gray-scale,  multiresohition  images  are  suitable  for  use  in  visual  amulators. 


INTRODUCTION  TO  WAVELETS 

Maiy  joint  position/spatial-fi-equency  techniques  have  been  devdoped  for  overcoming  the 
limitations  of  FA  described  above  (cf,  Jacobson  &  Wechslo’,  1988).  The  moa  straightforward 
way  to  obtain  a  joint  representation  is  to  simply  multiply  the  signal  by  an  appropriate  window 
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fiinctioii,  compute  the  Fourier  transform  of  ti»  product,  and  rqjeat  the  procedure  for  windows 
translated  in  position.  Tins  procedure  defims  the  short-time  (or  windowed)  Fotaier  transom 
that  allows  a  signal  (or  image)  to  be  rq)resented  in  both  po^on  and  qiatial  frequency.  Any 
smooth,  spatially-localized  function  can  be  used  as  a  window  hinction,  and  many  have  been  tried. 
Gabor  (1946)  has  shown,  howevo-,  that  a  window  fimction  in  the  form  of  a  gausaan  provides  the 
greatest  joint  localization  in  position  and  spatial  frequency.  The  Gabor  represoitadon  is  a  special 
case  of  what  are  known,  in  the  physics  lito’ature,  as  cohetont  stat»  associated  with  the 
Weyl-Heisoibeig  group  (Daubechies,  1992). 

In  the  presem  context,  the  We^-Heisoiba’g  group  represents  one  of  two  classes  of  what 
are  now  referred  to  as  wavelets.  Wavdets  are  sets  of  fimctions  that  are  formed  by  applying 
dilation  and/or  translation  operators.  Thus,  we  can  define  a  class  of  Weyl-Heisenberg  group 
wavdets,  vduch  are  generated  by  a  translation  in  both  ^ace  and  ^atial  frequency  as  follows; 


g»As)^g(s-b)  ef^  ,  (1) 

where  g(s-b)  is  a  window  function  that  can  be  trandated  in  space  by  an  amount  proportional  to  b, 
and  the  complex  exponential  represents  spectral  functions  translated  in  fi^uency  by  an  amount 
proportional  to  a.  Examples  of  Weyl-Hdsenberg  group  wavdets  are  shown  in  Figure  la.  These 
wavdets  all  have  the  same  window  function  (a  gaussian  in  this  case)  that  has  been  translated  in 
space  and  multiplied  by  spectral  functions  of  various  spatial  fivquendes. 

The  second  class  of  wavelets  are  the  c0rK  group  wavelets,  which  are  goierated  by 
dilations  and  translations  in  space  as  follows: 


v|/..D(s)=|ar5\|/(^^)  ,  (2) 

where  the  parameter  a  represents  dilation,  and  the  parameter  b  represents  translation.  Examples 
of  afiSne  group  wavelets  are  shown  in  Figure  lb.  These  wavelets  were  developed  to  analyze 
transients  in  signals  (or  images),  which  had  previously  been  dealt  with  by  non-optimal,  and  largely 
ad  hoc,  techniques. 
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(a)  Weyl-Heisenberg  Wavelets 


(b)  AfRne  Wavelets 


Figure  1 

Examples  of  (a)  Weyl-Heisenberg  Wavelets  and  (b)  AfBne  Wavelets. 
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As  noted  above,  one  nuyor  motivation  for  devdoping  wavdet  techniques  is  the 
requirement  for  a  joint  ^>ace  and  spatial  fiequency  analysis  to  overcome  the  limitation  of  FA, 
i^ch  assumes  infinite  resolution  in  one  domain  and  no  resolution  in  the  other.  In  tl^  context  of 
image  representation,  wavdets  are  used  to  decompose  images,  and  the  decomposition  is  most 
efficient  when  the  image  information  represented  by  each  of  the  wavelets  in  the  set  is 
nonredundant.  Wavdets  which  aUow  such  a  nonredtmdant  representation  are  called  orthogonal, 
and  can  be  used  to  represent  images  with  a  minimum  number  of  coefficients.  The  affine-group 
wavelets  are  orthogonal  whereas  the  Weyl-Hdsenberg  group  wavelets  are  not.  The  affine-group 
wavelets  are  also  better  suited  for  representing  transients  and  other  high-fi’equency  image 
componoits.  Since  our  miyor  objective  is  image  data-reduction  (as  a  first  step  toward  image 
compression),  we  are  more  concerned  with  the  above-mentioned  desirable  properties  of 
affine-group  wavdets  and  less  concerned  with  the  high  joint  localization  attainable  using  the 
Weyl-Hdsenb«g  group  wavelets.  We  will,  therefore,  limit  the  remainder  of  our  discussion  to 
affine  group  wavelets. 

CONSTRUCTION  OF  AFFINE-GROUP  WAVELETS 

The  affine-group  wavelets  to  be  described  here  fall  into  two  categories  that  will  be 
rrferred  to  as  scaling  functions  and  wavelets.  Scaling  fimctions  are  obtained  as  solutions  to  a 
two-scale  difference  equation  of  the  form: 

<Kx)  =  Scfc({)(2x-;0  ,  (3) 

k 


where  the  are  a  set  of  coefficients. 

Wavelets  are  obtained  as  solutions  to  a  difference  equation  of  the  form; 


V(x)  =  Z  H2x~k) 

k 


(4) 
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where  N  is  the  number  of  coefficients.  Unique  solutions  to  Equations  3  and  4  are  guaranteed  if 
tlw  sum  of  the  coefficients,  is  exactly  two  (cf ,  Stnmg,  1989).  The  solutions  can  be  obtained 
startup  from  any  of  a  number  of  initial  conditions,  denoted  for  ^x).  The  simplest  initial 
condition  is  the  "box"  function,  shown  in  Figure  2a.  Using  the  box  function,  various 

interesting  and  useful  scaling  functions  and  wavelet  pairs  can  be  obtained  by  using  different  sets  of 
coefficients.  For  instance,  when  Co=  c,  =  1,  a  scaling  function  can  be  obtained  as  a  solution  to 
Equation  3; 


♦‘i’W  =  +‘2’(2x)-Kj<?(2x-1)  ,  (5) 

which  is  simply  the  box  function  again  (i.e.,  This  construction  is  shown  gr^hically  in 

Figure  2.  The  assodated  wavelet,  for  Co=  c,  =  1,  is  the  solution  to  Equation  4; 

v‘i*(x)  =  -*?’(2x)  +  C(2x-1)  ,  (6) 

which  is  illustrated  graphically  in  Figure  2.  The  function,  yf/x),  known  as  Haar's  wavelet,  is 
simple  in  form  and  was  the  first  function  to  be  used  to  generate  an  orthogonal  basis.  However, 
whereas  y/f/x)  is  well  localized  in  space,  it  has  a  discontinuity,  and  so  does  not  provide  good 
frequency  localization.  It  is,  therefore,  not  ideal  for  image  representation. 

Another  scaling  function  and  wavelet  can  be  obtained  again  using  the  box  function,  but 
with  the  coeffidents  C(,=  Vi,  c,  =  1,  and  =  Vi.  The  scaling  function  is  obtained  using  Equation  3 
as: 


and  the  corresponding  wavelet  is  obtained  from  Eqn.  (4)  as: 


V?(x)  =  -^<2’(2x)  +  ((,<g'(2*-1)-l4.‘|’(2x-2)  .  (8) 
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original  box  I 

function  Haar  wavelet 


Figure  2 

Generation  of  a  Scaling  Function  and  Haar  Wavelet  (c)  from  the 
Scaled  and  Translated  "Box"  Functions  shown  in  (a)  and  (b). 
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The  resuhs  of  four  iterations  using  Equation  7  and  8  are  depicted  gr^hically  in  Figures  3a  and  3b, 
respectively.  The  scaling  ftinction  shown  in  Figure  3a  is  triangular  function,  while  the  wavelet  of 
Figure  3b  is  the  so-called  "hat"  function.  Althou^  the  hat  function  is  smoother  than  the  Haar 
wavelet,  it  is  not  orthogonal  to  its  translations  and  dilations  and  thus  yields  a  nonorthogonal  basis. 

A  final  example  of  a  useful  scaling  function  and  wavelet  pair  is  the  cubic  B-spline  obtained 
using  the  coefiBcients  Co=  1/8;  c,  =  4/8;  Cj=  6/8;  Cj=  4/8;  and  c^=  1/8.  The  first  through  fourth 
iterations  of  these  functions  are  shown  in  Figures  4a  and  4b,  respectiveh'.  The  scaling  function  is 
smooth  and  may  represent  a  good  approximation  to  a  gaussian  low-pass  filter.  The  wavelet,  on 
the  other  hand,  has  a  triphasic  form  which  is  similar  to  that  of  both  measured  and  theoretical 
visual  receptive  field  profiles  (cf ,  Young,  1987). 

A  Continuous,  Localized,  Orthonormal  Wavelet 

A  family  of  wavelets  (actually  scaling  functions  and  wavelets)  of  both  practical  and 
theoretical  importance  was  devised  by  Daubechies  (1988).  The  scaling  functions  and  wavelets 
shown  in  Figures  2  through  4  were  generated  using  rational  coefiBcients.  Daubechies  found  a  set 
of  irrational  coefiBcients,  which  resulted  in  scaling  functions  and  wavelets  that  are  orthonormal  to 
their  integer  translations  and  dilations.  These  functions  are  particularly  useful  because  both  they 
and  their  Fourier  transforms  have  compact  support  (i.e.,  are  of  finite  extent)  and  are  continuous. 
Since  we  will  be  using  this  scaling  function  and  wavelet  to  decompose  and  reconstruct  images,  the 
technique  for  generating  them  now  be  described  in  detail. 

Since  all  initial  conditions  converge  to  the  same  function,  for  simplicity  we  will  again  use 
the  box  function,  to  generate  the  scaling  function  and  wavelet.-  The  set  of  four  coefiBcioits 
suggested  by  Daubechies  (1988)  are; 

co  =  i(l  +  V3)  =0.683013 
Cl  =i(3  + 1/3)  =  1.18301 
C2  =  i(3 -1/3)  =0.316987 
C3  =  i(l-y3)  =-0.183013 

Using  Equation  3  and  following  the  procedures  described  earlier,  the  first-iteration  scaling 
function  is  obtained  as: 
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Scaling  Function  and  (b)  a  Cubic-Spline  Wavelet 


I 


(9) 


♦<i>(x)  =  0.683  ♦5>(2x)  +  1.183  ♦?(2x- 1) 

+  0.317  ^^(2x-2)  +  (-0.183)  ♦^(2x-3) 

The  function  is  shown  gi^hicatly  at  the  top  of  Figure  5a.  The  second-iteration  scaling 
function  is  then  obtained  from  as  follows: 


<|>?(x)  =  0.683  4><i>(2x)  +  1.183  ((>^i^(2x  - 1) 

+  0.317  (►§(2x-2)  +  (-0.183)  (►^J(2x-  3)  (10) 

This  procedure  is  repeated  to  obtain  all  higher-it«ation  scaling  fiincdons  as  requireo.  The  second 
dirough  fourth-iteration  Daubediies  scaling  functions  are  also  shown  in  Figure  Sa.  As  a  practical 
matto-,  no  significant  chaises  occur  in  this  &inction  after  the  seventh  ita'ation. 

In  analogous  ftstdon,  but  using  Equation  4  rather  than  Equation  3,  the  first-iteration 
Daubechies  wavelet  can  be  obtained  as  follows ; 

V  oV)  =  0.183  <(»?(2x)  +  0.317  ♦?(2x- 1) 

-1.183O2x-2)  +  0.683«j>^/(2x-3)  .  (11) 

The  first-iteration  wavelet  is  shown  in  Figure  Sb  along  with  the  second-  through  fourth-iteration 
wavelets.  In  the  next  section  we  will  describe  a  technique  for  constructing  an  orthonormal  basis 
from  the  scaling  function  and  wavelet  of  Figure  S.  Such  a  basis  yields  a  bett^  joint  representation 
than  does  that  associated  with  Haaris  wavelet,  but  at  the  expense  of  basis-function  regularity. 
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IMAGE  refresentahon  using  wavelets 


btrodvctioB 

Image  decompo^on  is  a  procedure  fi>r  obtaining  an  ahemative  representation  (i.e.,  a 
transformation)  of  an  image.  The  procedure  results  in  a  set  of  numbers,  called  expansion 
coefficietOs,  which  represoit  the  amjditudes  of  each  of  a  s^  of  q)ecialized  functions  (known 
coHectivdy  as  a  basis)  that  can  then  be  summated  to  produce  the  original  image.  The  usefolness 
of  image  decon^sition  is  that  certain  image  components,  i^ch  may  be  naore  important  than 
others  in  a  particular  context,  can  be  idoitified  and  isolated.  Once  isolated,  these  components  can 
be  selectivdy  subjected  to  further  procesting.  In  the  cmitext  of  inoages,  the  decomposition 
procedure  partitions  information  into  two  domains,  namdiy  space  (or  portion)  and 
spatud-frequency:  Since  the  space  domain  is  the  image  itself  most  of  the  novel  information 
resulting  from  a  decomposition  is  in  the  spatial-frequoicy  dcnnain.  Further,  it  is  often  the  case 
that  componoits  of  interest  are  segregated  in  the  ^atial-frequency  domain.  This  observation  is 
the  m^or  motivation  for  the  development  of  multiresolutton  (or  pyramid)  decompotition 
techniques. 

In  multiresohition  decompotition  (cf ,  Akansu  &  Haddad,  1992),  an  image  is  filtered  into 
low-  and  high-frequency  components  and  subsampled,  and  this  procedure  is  repeated  as  required. 
Multiresolution  analytis  was  originally  done  (Burt  &  Adelson,  1983)  using  second  derivatives  of 
gausaans  (i.e.,  Liq>lacian  functions).  Mallat  (1989)  formalized  the  procedure  and  extended  it  by 
using  symmetric,  non-orthogonal  wavelets,  since  orthonormal  wavdets  were  not  available  at  that 
time  (see  Figure  1  in  his  ^pendix  A).  As  noted  above,  Daubechies  (1988)  has  derived 
orthonormal  wavelets,  which  while  not  symmetnc,  do  have  compact  support  and  reasonably  good 
joint  localization  in  space  and  spatial-frequency.  It  was  not  immediately  clear  that  orthonormal 
wavelets  were  suitable  for  multiresolution  image  decompotition  and  reconstruction  because  they 
were  not  symmetric,  and  image  processing  was  usually  performed  with  symmetric  kernels. 
However,  it  is  now  well  established  that  Daubechies*  wavelets  are  iq>propriate  for  this  purpose 
(Akansu  &  Haddad,  1992).  Further,  they  are  computationally  eflScient  in  that  decomposition  and 
reconstruction  can  be  done  using  filters  constructed  from  as  few  as  four  coefficients.  In  the  next 
section  we  will  outline  a  method  for  decompoting  and  reconstructing  an  image  using  the 
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cooqxrtatioiiid  {NiDcedures  desoflMd  by  MaB«t  (1989),  and  the  wavdet  basts  derived  by 
Daubechies  (1988). 


Ibu^  DecompositioB 

In  order  to  p^orm  a  imiltiresolution  image  deconqtosition  using  the  scaling  functions  and 
wavelets  described  eaiiio*,  we  use  ^j/x)  and  y/p(x)  (see  Figure  S),  and  define  the  following 
orthoirior7fta/>wive/erdas»(Daubechie^  1988): 


<|>>,(x)  =  2^(|iD(2^x-n) 

(12) 

M/>,(x)  =  2-^\l/o(2^x-n)  . 

(13) 

We  can  fiirtho'  define  (cf.,  Daubechies,  1988;  Mallat,  1989)  filters  h  and  g  as  the  inner  product  of 
particular  scaling  functions  and  wavelets  as  follows: 

/i(n-2/^  =  2-'«J*(|)(Kx-(n-2IO)<fa,  (14) 

and 

gin-2k)  =  2-^^jyvii)¥x-{n-2l()dx.  (15) 


The  functions  h(n-2k)  and  g(n-2k)  can  be  conadered  the  pulse  responses  of  low-pass  and 
high-pass  filters,  respectively.  Thus,  using  the  multiresolution  technique  developed  by  Mallat 
(1989)  for  symmetrical  wavelets,  we  can  expand  a  one-dimensional  signal,  S  (here  denoted  5^, 
into  two  components,  and  S^,  as  follows: 


Si.1  =  h(n~2k)  •  Su> 

ft 

and 

Shi  =  l,g(n-2k)SLo 

ft 


(16) 

(17) 
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Thin,  Su  and  Sg,,  i^ch  are  the  low-pass  and  U^pass  versions,  respectiveiy,  of  the  image, 
are  obtained  by  convolving  .S'^  with  the  low-pass  and  hi^pass  filters  h(n)  and  g(n).  The  fiu:tor 
of  two  in  the  argumoit  of  both  Equations  16  and  17  indicates  that  the  filter  (conv(duti<Mi)  output 
is  subsanqried  by  a  fiictOT  of  two. 

The  deconqwsition  {nocess  of  Erputions  16  and  17  can  be  rqjresented  as: 

S=Sija->SLt,Sm  , 

indicating  that  5*^  is  decomposed  into  S^,  and  where  represents  a  lower  resolution 
(smoother)  verskm  of  and  Sg,  rqnesoits  the  details  of  that  are  misang  fix>m  S^.  The 
decomposition  process  can  be  continued  by  finther  decmnposing  such  that: 


Sli  ->  St2  ,  Sh2  , 


and  more  generally. 


Sip  — ►  Sip^i  ,  (18) 

^^wre,  again,  is  the  low-pass  filtered  version  of  5^  and  Sg^j  is  the  high-pass  fiHaed 
veraon  of  iS'^.  Note  that  if  h(n)  is  designed  to  be  symmetric  [i.e.,  h(n)  =  h(-n)]y  then  Equations  16 
and  17  can  be  written  as: 


and 


n 


S»^i(K)-£s('l)-Si,(n-2A)  , 

« 


(19) 


(20) 
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req)ectivdy.  The  operatkuu  of  Elutions  19  and  20  are  equivaloit  to  taldng  every  other  sano^le 
fixnn  the  result  of  crmvdving  5^  with  A  and^. 

In  sununaiy,  the  decon^IXMition  of  a  sigoid,  V.  into  Hs  con^Kments,  5^  and  is 
adoeved  by  cotivolviiig  that  signal  with  the  filters  h  and  g,  rejq)ectivdy,  fi)llowed  by  dedmating 
the  resultant  signals  by  a  fiu^tor  of  two  0.e.,  taldng  every  other  san^de).  A  kvd-one 
<toc(mq;)osition  of  the  signal  5^,,  ndiose  lerigth  is  AT,  results  in  the  components,  5^  and  eadiof 
lex^giikN/2.  Thif  decomposition  process  is  depicted  in  Figure  6.  Unearty  crmvolving  ndiidi 
is  of  length  with  dther  h  or  g;  idiidi  are  of  length  L,  results  in  a  signal  of  length  N-^L-I.  To 
deal  with  this  problem  we  used  circular  convolution.  As  Equations  16  and  17  indicate,  if  the 
signal  5*^  is  of  length  then  its  comptments,  and  will  each  be  of  length  N/2*. 

Shown  in  I^gure  7  is  a  sdiematic  for  extending  the  deconq>osition  process  to 
two-dimensional  signals  such  as  images.  Figure  8  gives  a  more  detailed  graphical  description  of 
the  two-dimensional  decomposition  process.  Signals  correq^nding  to  and  S^,  are  drtained  by 
filtering  eadi  row  of  the  original  image,  and  are  denoted  LP  and  HP  in  Figure  8.  These  signals 
form  the  intermediate  image.  The  filterirg  p/ocess  is  ttoi  rq)eated  for  the  columns  of  the 
intermediate  image  therdiy  producing  the  levd-one  decomposition.  An  example  of  the  levd-mie 
decomposition  applied  to  a  real  image  is  shown  in  I^giure  9.  In  this  figure  (as  wdl  as  in  Hgure  10 
below),  the  pixd  values  in  each  of  the  three  cdls  containing  high-pass  (HP)  information  have  been 
artificial^  increased  to  make  than  more  easily  visible.  An  exanople  of  a  hi^ier-order  (level-three) 
decompodtion  is  shown  in  Figure  10. 

Image  RecoBStmetioB 

In  the  presoit  context,  image  reconstruction  is  the  reverse  operation  of  image 
deconqx>sition.  The  higher  resolution  agnal,  can  be  reconstructed  fium  its  components, 
and  as  follows: 


St^i  =  £  h{n-K)  SL  +  Z 

K  k 


=  S'Lp*h(n)  +  Si^*g(ri)  , 


(21) 


15 


16 


I 


< 
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Figure  9 

An  Example  of  a  Level-One  Decomposition 
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whov  (♦)  is  the  convolution  operator  and  is  an  interpolated  version  of  (i.e.,  with  a  mean 
luminance  value  inserted  after  every  sample  of  Schematics  for  the  one-  and  two-dimoisional 
reconstruction  process  are  shown  in  Figures  1 1  and  12,  respectively. 

ImplemeiiUitioii 

It  follows  from  Daubechies  (1988)  that  g(n)  and  h(n)  are  related  as; 

g(n)  =  (-1)«/i(1-n) 

vdiere  h(J-n)  is  the  mirror  filter  of  h(n).  Thus,  as  a  practical  matter,  if  g(n)  is  known,  h(n)  can  be 
obtained  directly.  We  will  now  show,  further,  that  h(n)  and  g(n)  are  equivalent  to  the  c^'s  of 
Equations  3  and  4.  The  significance  of  this  equivalence  is  that  it  allows  the  required  expansion 
coefficients  to  be  computed  from  the  scaling-function  coefficients,  c^.  Thus,  it  is  not  necessary  to 
generate  the  full  wavelet  waveform  or  to  use  that  waveform  in  the  computation  of  the  inner 
products  from  which  the  expansion  coefficients  are  obtained. 

As  implied  by  Equations  16,  17,  and  18,  we  can  decompose  an  image,  H,  at  resolution 
level  T'  with  respect  to  the  functions  and  y.  Those  equations  can  be  rewritten  as: 

Sip  =  H(x)athv9l2^^  =  Z  .<!>)«>  +  E  ,  \|//k>  (22) 

k  k 

where  and  v|/j,,  are  the  low-  and  high-pass  filter  coefficients  (or  basis  functions),  and  so  the 
summations  represent  low-  and  high-pass  filtering.  The  inner  products  of  Equation  22  give  the 
wavelet  (or  decomposition)  coefficients.  For  instance,  when  extended  to  two-dimensions,  the 
first  inner  product  is  simply  the  upper  left  comer  of  Figure  9.  We  will  now  show  that  the  inner 
products  in  the  summations  are  equivalent  to  the  c^'s  of  Equation  1. 

We  begin  by  using  Equations  4  and  13.  Specifically,  from  Equation  4  we  have; 
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23 


Figure  12 

Schmatic  of  the  Two-Dimensional  Reconstruction  Process 


y^-<x  -  n)  =  2  (-1)"c»^hi,  ♦C  •  C-'x  -  n)  -  *) 

k 

=  2  .  X  -  (2/1  +  A))  (23) 

* 

Substhutii^  Equation  23  into  Equaticni  13  ^es: 

yM*)=2(-1)“c*<».i)  •  2-4  •♦(2-»-’>x-(2n+(()) 

Jr 

=  2(-1)**'c»<h„-2-4-2'^+(2^^'>x-(2n+*))  , 

k 

and  usii%  the  notation  of  Equation  12  we  get: 

=  S  (-1)*^^  CN^k¥\)  •  2'^  ♦(^iwarHioOf)  (24) 

k 

Denoting  /  =  2w+A  and  k+I  =  (l-2n+l)  gives: 


VjnC^)  =  L  (-1)^^^Ca^2»H)  '  2"2 


and  denoting. 

g(J-2ri)  =  •  2"^ 

(25) 

^es. 

yVjn(x)  =  Xg(!-2ri)-^(f-i)j  . 

(26) 

The  high-pass  wavdet  coefiBdents  at  a  given  resolution  level  can  be  obtained  by  taking  the  inno- 
product  of  the  above  equation  with  the  image,  H : 
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(27) 


# 

<W ,  vjb>  =  2  2n)  (H , 

^lere  g  is  a  high-pass  fiher.  The  inner  prochict  (<  ~  >)  on  the  right  side  of  Equation  27 
represems  the  wavdet  co^dents  at  one  resolution  levd,  whereas  the  inna*  product  on  the  Idi 
side  of  the  equation  rq>resents  the  wavdet  coeffidents  at  the  next  highor  levd. 

Following  the  same  procedure  as  indicated  in  Equations  23  through  26,  the  low-pass 
wavdets  coeffidoits  can  be  obtained  as; 

.  <H>>,)  =  2/K/-2n).<H.*cMy>  .  (28) 

I 

where,  in  analogy  with  Equation  25,  hi  =  Citj2  . 

Equations  14,  IS,  and  22  diow  how  to  compute  a  coarser  approximatimi  from  a  finer  (me, 
as  wdl  as  how  to  compute  the  difiGvence  in  informati<m  betweoi  the  successive  approximations. 
To  (4>tain  a  finer  re|»esaitation  from  a  coarser  one,  we  proceed  sinularfy  obtaining: 

=  +  v/r>  (29) 

/  / 


SUMMARY 

The  essence  of  multiresoiution  image  decomposition  and  reconstruction  is  the  use  of  a 
pair  of  filterS"One  low-pass  and  one  high-pass.  The  filters  are  constructed  from  orthonormal 
wavelets  by  taking  the  inner  product  of  the  dilated  wavelet  and  a  shifted  version  of  the  undilated 
wavelet.  Image  decomposition  is  performed  by  fii^  pasdng  each  image  row  and  coiunm  through 
both  a  low-pass  filter  and  a  high-pass  filter.  Since  this  filtering  operation  reduces  the  image 
bandwidth,  fewer  pixels  can  now  be  used  to  represent  each  line  with  no  loss  of  information. 
Therefore,  the  output  from  each  filter  is  subsampled  by  eliminating  every  other  pixel.  The 
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|HX)cedure  is  rq)eated  for  each  row  and  column  produdng  four  images  each  one^uarter  the  size 
of  the  original.  The  procedure  is  q>plied  recursivdy  to  obtain  additional  levds  of  decomposition. 

Deconq)osition  is,  as  a  practical  matter,  always  performed  down  to  the  2  x  2  pixel  levd. 
After  deconqmation,  redundam  information  can  be  removed  up  to  any  chosen  level. 
Reconstructions  can  thm  be  pofonned  up  to  various  levels  in  order  to  determine  the  lowest  level 
^di  gives  an  accq>table  reconstruction.  We  start  the  reconstruction  from  the  lowest 
decomposition  level  (the  2x2  level).  The  subimage  corresponding  to  that  level  is  first  processed 
through  the  same  high'pass  and  low-pass  fibers  used  in  the  decomposition.  The  resulting  filtered 
images  are  then  doubled  in  size  by  inserting  the  mean  luminance  value  between  all  pixels,  and 
corresponding  pixel  values  are  added  together  to  produce  a  single  inuige.  This  procedure  is  then 
repeated  on  this  image  and  all  subsequent  images  until  the  original  image  size  is  achieved.  The 
original  image  can  be  reconstructed  exactly  by  this  procedure.  The  point  of  the  multiresolution 
analysis,  in  the  present  context,  is  to  identify  information  that  can  be  r«noved  without 
significantly  affecting  the  appearance  of  the  reconstructed  image.  The  advantage  of  the 
orthonormal  wavelet  technique  is  that  it  provides  a  formal  method  for  reconstructing  an  original 
image  from  the  multiresolution  decompotition.  Prior  to  the  development  of  this  tedinique, 
multiresolution  reconstruction  could  not  be  performed  error  fi^  and  not  so  simply  as  in  the 
orthonormal  wavelet  case  wherein  reconstruction  is  simply  the  reverse  process  of  decomposition. 
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APPENDIX 


Source  code  for  the  fvogram  wavelts.c  v^ddi  performs  image  deconqwsHioii  and 
leconstniction  with  respect  to  a  contiiiuous,  orthonormal  wavdet  basis. 
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